How small can thermal machines be? Towards the smallest possible refrigerator 
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We investigate the fundamental dimensional limits to themodynamic machines. In particular we show that 
it is possible to construct self-contained refrigerators (i.e. not requiring external sources of work) consisting of 
only a small number of qubits and/or qutrits. We present three different models, one consisting of two qubits, 
one consisting of a qubit and a qutrit with nearest-neighbour interactions, and one consisting of a single qutrit. 
We then investigate fundamental limits to their performance; in particular we show that it is possible to cool 
towards absolute zero. 



When over two centuries ago Sadi Carnot |1| set out to 
study the physics of steam engines - and in the process es- 
tablished thermodynamics, one of the cornerstones of physics 
- the key to his progress was to abstract from real machines 
to idealised, "model independent" heat machines. What he 
found is that although the properties of each machine depend 
on the details of its construction, the fundamental limit to their 
efficiency is independent on these details. But can physics in- 
deed be completely left out? Here we return to physics and 
ask about other fundamental limits, specifically, how small 
can thermal machines be? Is there a fundamental limit to their 
size? And, when they are very small, are there additional con- 
straints on their performance? Can small machines work as 
well as macroscopic machines, or is there a complementarity 
between size and performance? For example, can small ma- 
chines be constructed that can cool arbitrarily close to absolute 
zero, or does size impose a fundamental limit? 

In the present paper we approach these questions in the 
framework of quantum mechanics. Most importantly, quan- 
tum mechanics provides a natural and universal notion of 
"size", namely the dimension of the Hilbert space that de- 
scribes the system. It is this measure of size that we will use 
here, to characterise the thermal machines, not their extent in 
space or their mass. Nevertheless, it is very often the case that 
small Hilbert space dimensions are associated to microscopic 
systems, (atoms, spins, etc.) and hence with small space di- 
mensions and small mass. 

The particular thermodynamic machines we study here are 
"self-contained" refrigerators. By "self-contained" we mean 
that they do not use any outside source of work. That cooling 
can be achieved when work is performed on a microscopic 
system is well-known. Indeed, laser cooling H IH is the 
work-horse of the entire field of atomic and molecular op- 
tics. Such externally driven refrigerators have also been stud- 
ied theoretically, where they have received much attention re- 
cently, e.g. i5ji6i[7iia[3iTni[iii[i2[i3i[ii[i5i[ig 

However, in this instance one can hardly call the microscopic 
system that is getting cooled a self-contained fridge - one also 
needs, for example, a macroscopic laser. In this paper we are 
interesting in accounting for all the degrees of freedom that 
constitute our fridge; this is the reason for focussing on ma- 
chines that do not require external work. On the other hand, of 
course, no thermal engine can work without a source of free 
energy. As in all thermal engines we provide this by allowing 



the fridge to be in contact with two heat reservoirs at different 
temperatures Tc < Th. 

Apart from the implications for fundamental physics, we 
envisage that this work may be relevant to a number of fields. 
One such field is biology. Here, for example, if cooling of the 
active site of a protein can be achieved then increased catalysis 
rates may be possible |18|. It is an intriguing question as to 
whether simple mechanisms, such as those we describe in this 
paper, are used by biological systems. A second field is the 
rapidly developing field of nanotechnology where the benefits 
from cooling at the atomic scale are clear. 

In this work we will present three models. The first re- 
frigerator consist of two qubits (two level systems). Building 
upon this idea we present a second fridge, consisting of a qubit 
and a qutrit (three level system) and where all interactions are 
nearest neighbour. As a last example we present a refrigerator 
consisting of a single qutrit arguably the smallest refrigerator 
possible. Finally we prove that there are no fundamental lim- 
its to how close to absolute zero such small refrigerators can 
cool. 



1. MODEL I: TWO QUBITS 

The first model we present consists of three qubits, two con- 
stituting the refrigerator, and one being the object to be cooled. 
It is inspired by the method of algorithmic cooling [19J, par- 
ticularly by the few qubit version in 1201 . 



A. Functioning principle 

Consider first two qubits, immersed in a bath at temperature 
Tc (the 'cold' bath). Qubit 1 is the object to be cooled, while 
qubit 2 will eventually play the role of the "spiral" that takes 
heat from qubit 1 and dissipates it into the environment. 

The "free" Hamiltonian Hq for the two qubits is 

Ho = ^iH^^) +^2n(^) (1) 

where H*^*^ = is a projector for qubit i, and we denote 

by |0) and |1) the ground state and excited states respectively. 
The energies of the excited states of the two qubits are Ei and 
E2 respectively, while their ground state energies are zero. 
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We take E2 > Ei, hence qubit 2 has a larger energy level 
separation than qubit 1 . 

At equilibrium each qubit is in a thermal state, r^, which 
has the form 

= iVie-^^n'''/'=^= (2) 

where A^^ = (1 + 6"^*/^^^=) ^ and A: is Boltzmann's constant. 
Since the qubits do not interact the total thermal state is simply 
the direct product state 

P12 Ti r2 

A convenient way to represent the thermal state ^ is in terms 
of the probabilities, and 1 — r^, to find the qubit in the 
eigenstates |0) and |1) respectively, 

r^ = r,|0)(0| + (l-r,)|l)(l|, (4) 

where = 1/A^i. 

We will denote by Tf^ the equilibrium temperature of each 
qubit when refrigeration is occurring. Cooling of qubit 1 
means reaching a temperature T^^ < Tc, which corresponds 
to a larger probability rf ^ > ri of being in the ground state. 

The idea of algorithmic cooling is to increase ri by trans- 
ferring excitations to the second qubit. Specifically, two par- 
ticular eigenstates of ([s]) are 1 10) and |01). The former appears 
with probability (1 — ri)r2 while the latter with ri(l — r2). 
Suppose now that we perform a unitary SWAP operation U, 
which swaps these two states, 

|10) ^ |01) (5) 
|01) ^ |10) (6) 

while leaving the others unchanged. Since we take E2 > Ei, 
the probability (1 — ri)r2 is larger than (1 — r2)ri, that is 
we are more likely to be in the state |10) (where qubit 1 is 
excited) than in |01) (where it is in the ground state). Thus 
after the swap we increase the probability of qubit 1 being in 
the ground state (we cool it) and we increase the probability 
of qubit 2 being in the excited state (we heat it). 

Since the qubits are in contact with a thermal bath, if we 
don't do anything else, the qubits will eventually return to 
the environmental temperature Tc. However, if we keep on 
repeating the unitary transformation U, qubit 1 will reach an 
equilibrium temperature T^^ lower than Tc and qubit 2 a tem- 
perature larger than Tc. 

The procedure described above however requires external 
work to be performed on the qubits. Indeed, the energy E2, 
of 1 01), is larger than the energy Ei, of |10). This is accom- 
plished by some external system that induces the unitary trans- 
formation U ; for example, in an NMR experiment this could 
be done via a sequence of pulses of an external magnetic field. 

The idea behind our fridge model is then to replace the ex- 
ternal work by a different source of free energy. Specifically, 
free energy can be obtained whenever a system has access to 
two thermal baths at different temperatures. In our model this 



is accomplished by adding a third qubit, in contact with a ther- 
mal bath (the 'hot' bath) at a higher temperature > Tc. 
This qubit plays the role of the "engine". Thus our fridge 
consists in qubits 2 and 3 (the spiral and engine respectively); 
qubit 1 is the object to be cooled. 

To enable transitions between different states without an in- 
put of external energy, we take the engine qubit to have the 
energy level spacing Es, such that Es = E2 — Ei. With 
this condition we now have two degenerate energy eigenstates 
1 010) and 1 101) and hence we can swap the two states without 
requiring work. The interchange 

|101) ^ |010) (7) 
|010) ^ |101) (8) 

accomplishes on qubits 1 and 2 the transformations ([5j and 

For transitions between |010) and |101) to occur we intro- 
duce an interaction Hamiltonian 

^znt=^(|010)(101| + |101)(010|) (9) 

where g determines the strength of the interaction. 

Furthermore, for simplicity, we will assume that this inter- 
action is taken to be weak compared to the energy level spac- 
ing, Ei ^ g. In this regime the interaction will not signifi- 
cantly alter the energy eigenvalues or eigenvectors of the sys- 
tem (which remains governed by Hq) and hence we can mean- 
ingfully talk about the temperature of the individual qubits, as 
each qubit will remain in the standard form of a thermal state 
([2]), with El and E2 the same as they were in the absence of 
interaction. 

To understand the functioning of the fridge first note that 
the interaction Hamiltonian can swap, without impediment, 
the states |010) and |101), since they are degenerate in energy 
in-so-far as the free Hamiltonian is concerned. However, if all 
qubits were kept at the same temperature, the system would be 
at equilibrium since the the probability of the flip ^ is equal 
to that of ([8]). In order to drive the transitions in one direction, 
such as to cool qubit 1, we place the third qubit in contact 
with a hotter bath. In this way the probability of the state 
1 101) becomes larger than that of |010) and so we enhance 
the probability of the flip ^ and diminish that of ([8]). It is this 
biasing of the interaction which has the effect of taking heat 
from qubit 1 into qubit 2, creating a refrigerator. 

In the next section we will describe our model in detail. 



B. Master Equation 

In the absence of any interaction, the three qubits have the 
free Hamiltonian 

Ho = EiU^^^ + E2U^'^^ + EsU^^^ (10) 

The interaction between the three particles, as previous stated, 
will be generated by the interaction Hamiltonian ^ 

We further need to model the fact that each qubit is taken 
to be in contact with a thermal bath. The simplest model of 
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this scenario is to imagine that with some probabiHty density 
Pi per unit time that each qubit may be thermaUsed back to 
its initial thermal state ([2]). (Note that for qubit 3 the bath 
temperature is Th and not Tc). Mathematically we model this 
by the non-unitary evolution 



(11) 



Here Pi quantifies how well insulated each particle is relative 
to the bath; when pi is small the particle is hardly ever ther- 
malised by the bath. 

All together this leads to the master equation 



dp 
dt 



-i[Ho ^ Hint, p] ^^PiinTiip- p), (12) 



(see Appendix A). Note that ([12]) can easily be rewritten ex- 
plicitly in Lindblad form. 



C. Equilibrium Solution 

We will be interested in solving for the equilibrium state, 
PEq, which satisfies 

3 

= -i[Ho + Hint, pEq] + ^PiiriTliPEq - pEq)- (13) 

i=l 

The equilibrium solution can be easily found analytically, 
but it has a complicated dependence on all the parameters, so 
it is much more illuminating to present a numerical analysis. 

Figure [T] shows the dependence of the temperature of 
qubit 1 (the system to be cooled) as a function of the temper- 
ature of the hot bath, T^. We see that when the temperature 
of the hot bath is larger than that of the cold bath, > 
i.e. in the normal regime of functioning, qubit 1 is cooled to a 
temperature T^^ < hence our refrigerator indeed cools. 




FIG. 1: Equilibrium temperature of qubit 1, Tf^, in contact with 
a bath at temperature Tc = 1, versus bath temperature Th of qubit 
3, with all insulation parameters taken to be equal, p^ = p. The 
different curves correspond to different values of g/p. We take the 
energy level spacings Ei = E3 = 1,^2 = 2. Inset: Schematic 
diagram showing energy levels and interaction. 



We notice that when = that is, when we supply no 
free energy, the temperature of qubit 1 is unchanged. Finally 
we also notice that if the temperature of the "hot" bath, is 
taken to be actually smaller than that of the "cold" bath, Tc, 
than our refrigerator works effectively in reverse and it warms 
qubit 1 to a temperature higher than the hottest environmental 
temperature (which in this case is Tc). 



D. Parameter dependence and Zeno effects 

A natural question to ask is how the behaviour of the fridge 
changes as we vary the parameters pi independently. Indeed, 
in a 'standard' refrigerator we do not want all parts to interact 
with the environment equally: the inside of the fridge has to 
be well insulated to maintain a low temperature while the spi- 
ral at the back of the fridge has to interact strongly with the 
environment so that it dissipates heat quickly. 

Concentrating first on qubit 2, (i.e. the spiral) we expect 
that as p2 becomes larger the performance of the fridge should 
increase (T^^ should decrease). Furthermore, qubit 3 plays 
the role of the "engine" of the refrigerator, which it achieves 
by pumping heat from the hot environment into the system. 
We expect that the best performance is therefore achieved 
when it interacts strongly with its environment, as this allows 
it to extract heat at the highest rate. 

In Fig. [2] we plot separately the dependence of T^^^ on p2 
and ps. For small values of p2 and ps we indeed observe the 
expected behaviour, however as we increase them further the 
performance degrades. The reason behind this is that there are 
quantum effects which come into play in the regime of strong 
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FIG. 2: Equilibrium temperature T^^ of qubit 1 as a function of 
the insulation parameter: (a) p2 and (b) ps. When p2 or ps vanish 
then we are unable to achieve cooling. For large values of p2 or ps 
the performance of the fridge degrade - in this regime a Zeno effect 
inhibits the functioning of the refrigerator. 
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coupling between the qubits and the environment. 

Specifically, the effect that we are observing is an instance 
of a Zeno effect II2TII . Thermalising each qubit is closely re- 
lated to a measurement; it is as though the environment 'mea- 
sures' the qubit as it puts it in a thermal state. Now, as we 
increase the rate of thermalisation we soon enter a regime 
whereby the interaction Hamiltonian Hint will not have time 
to work between successive thermalisations and hence the 
whole refrigerator is no longer able to function. 



E. Approaching absolute zero 

We would now like to look more closely at the performance 
of our refrigerator. A first important question is whether or not 
there slyq fundamental limitations to the temperature to which 
we are able to cool down the first qubit. We show that there is 
no such limitation. 

The minimal achievable temperature is limited by two ef- 
fects: a heat current flowing into the fridge due to imperfect 
insulation, and the actual cooling ability (i.e. the ability to 
cool given perfect insulation). It is the second aspect which 
we are interested in. 

We fix El, since this is a characteristic of the object to be 
cooled and not of the refrigerator and also fix Tc, as this is 
the environmental temperature. We increase E2 and such 
that the ratio Es/T^ remains constant and much less than 1. 
This results in increasing the population of the ground state 
for qubit 2, while maintaining a high probability of finding 
qubit 3 in the excited state. Put together these two elements 
mean that the interaction ^ becomes ever more biased as we 
increase E2. This leads to cooling as close as we want to 
absolute zero, as seen in Fig. [3] 

2. MODEL II: ONE QUBIT, ONE QUTRIT 

One drawback of the previous model is that the interaction 
Hamiltonian ([9]) is a three body interaction. Here we present 
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FIG. 3: Equilibrium temperature Tf^ of qubit 1 as a function of the 
energy level spacing E2 . As we increase the spacing we are able to 
cool the spin arbitrarily close to absolute zero. 



a model with only two-body nearest neighbour interactions. 

The model consists of three particles, where one is to be 
cooled and two to construct the fridge. Particle 1 and 3 are 
qubits and particle 2 is a qutrit (three level system). The en- 
ergy levels of each particle are as depicted in Fig. |4] (inset) and 
are such that the energy eigenstates |020) and 1 101) are degen- 
erate in energy. By introducing an interaction which can take 
the population of the latter into the former we can cool down 
qubit 1. 

We do this by introducing two separate interactions be- 
tween the particles via the Hamiltonians 

= 5(|02)(11| + |11)(02|)®]1(3) (14) 

= MW®(|01){10| + |10)(01|) (15) 

Neither ([14]) or ( p3] ) induces transitions between the two de- 
sired states. However, ([14]) causes transitions between |020) 
and 1 110) and ( [T5] ) between |110) and |101). Therefore, in 
second order we indeed induce the desired transition. 

Finally we bias this interaction as in the previous model, by 
taking particles 1 and 2 to be in contact with a bath at temper- 
ature Tc and qubit 3 at T^. 

The details of this model are given in Appendix B. Here 
we simply display the results obtained in Fig. |4] We observe 
that this refrigerator behaves qualitatively the same as in our 
previous model. 
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FIG. 4: Equilibrium temperature ^ of qubit 1 as a function of the 
bath temperature Th in which qubit 3 is immersed. We have taken 
Pi — P2 — P2, — Tc — 1 and E — 1. Inset: Schematic diagram 
showing energy levels and interaction. 



3. THE SMALLEST POSSIBLE FRIDGE 

In the previous model particle 3 is required only so that we 
can heat the lower transition on qubit 2 (between |0) and |1)) 
independently of the upper transition (between |1) and |2)) 
as we have taken the thermalisation to thermalise both transi- 
tions together. However, we could conceive a situation where 
each transition is thermalised independently, e.g. if each en- 
ergy eigenstate has a different spatial distribution then each 
one can be in contact with an environment at differing tem- 
peratures. In this situation it is then possible to construct a 
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FIG. 5: Schematic diagram showing a fridge consisting of a single 
qutrit (particle 2) along with the object to be cooled (particle 1) 



smaller refrigerator, by discarding the third qubit - the fridge 
now contains only a single qutrit and the interaction ([14]) (See 
Fig.|5]). We believe this is the smallest possible system which 
may be called a refrigerator. 



4. CONCLUSIONS 

We have presented three simple models which demonstrate 
that there is no fundamental difficulty in constructing small, 
self contained systems which act as refrigerators. We re- 
call that by 'small' we mean precisely that we take quantum 
systems consisting of very few states. Furthermore we have 
shown that it is possible to cool towards absolute zero. 

There are however many interesting questions and direc- 
tions which we have not yet even begun to explore, of which 
we shall mention a few. The first is to explore what can be 
said about the efficiency of the refrigerators we presented. Is 
it possible to quantify the efficiency? If so is our construction 



the most efficient or are there other Hamiltonian's which are 
better for cooling? 

Moreover, it is fundamental to ask whether or not there 
exists a complementarity between small dimension and effi- 
ciency - is it the case that you can only be large and efficient 
or small and inefficient? In other words, can a small machine 
reach the efficiency of an ideal Camot engine? Our particu- 
lar models do not reach this efficiency; indeed, both the spi- 
ral and the engine qubits reach equilibrium temperatures that 
are different by a finite (instead of infinitesimal) amount from 
the temperatures of their environments, which leads to irre- 
versible heat exchanges. However, could it be possible to find 
a better model? This is not clear, since in a Carnot cycle the 
system transitions through very many states. Here, however, 
it is not the case that we can pass through many states in a 
single cycle of refrigeration. The question is whether or not 
this affects the achievable efficiency of the refrigerator. 

Here we focused on refrigerators, but studying other ther- 
modynamic machines would be an interesting direction to ex- 
plore. Given a supply of free energy we have already shown 
how to cool a system, but are there other interesting tasks that 
we can also achieve given only an interaction with two ther- 
mal baths? Could we somehow extract the 'work' and use it 
to drive other molecular machines? 
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APPENDIX A: DERIVING MASTER EQUATION 

To find the master equation consider a small time interval 
5t around a time t^. To first order in 5t, the evolution of the 
density matrix p is given by 

p{tQ^5t) = (l-(5t(pi+p2+P3))p(^o) 

+ (^t(piriTrip(to)+P2r2Tr2p(to) (Al) 
+ P3T3Tr3p(to)) -i3t[HQ^ Hint,p{tQ)\. 

from which it follows that 

-£ = -i[Ho + Hint, p] + ^PiinTiip - p), (A2) 
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Note that ^V2) can easily be rewritten explicitly in Lindblad 
form. 



APPENDIX B: DETAILS FOR TWO-QUBIT-ONE-QUTRIT 
REFRIGERATOR 

The free Hamiltonian of the 3 qubits is taken to be 

Ho = EiU^^^ + ^2nf ^ + (^1 + E2)U^2^ + ^2nf ^ (Bl) 

where u['^ = and n^'^ = \2)i{2\. The interaction 

Hamiltonian is now 

H,nt = ^(|O2)(ll| + |ll)(O2|)0]l(^) 

+ /il^^) (|01)(10| + |10)(01|) (B2) 

where g and h are the coupling constants for each interac- 
tion. We again assume that the interaction strength between 
the qubits is small in comparison to the energy level spacing, 
Ei, that is g^h <C Ei. Each particle again interacts with a 
thermal bath. For qubits 1 and 3, the thermal states are given, 
as previously, by For particle 2 the thermal state is now 
given by 

r2 = exp(-(£;2nf ^ + {E, + E2)U^2^)/kT,) (B3) 

where = (1 + e"^^/^^^ + e-^E,+^,)/kT^yi ^ 
The master equation governing the evolution is 

= -i[Ho + H,nt. p] + ^P^(r,Tr,p - p) (B4) 

i=l 

for which we will be interested in solving for the equilib- 
rium solution pEq- As before, we will not present an analytic 
form for pEq but only study properties of it numerically. Fig- 
ure]?] displays the dependence of the equilibrium temperature 
of qubit 1, T^^ on T/^, the temperature of the hot bath, for 
various values of the parameters g/p and h/p. 



